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1. Introduction 

Entanglement is a quantum correlation that appears in composite systems and 
constitutes one of the main resources in quantum information [1]. In optics, parametric 
down conversion (PDC) in a non-hnear crystal [21 [3l H] is the standard technique to 
generate entangled photon pairs. However, the statistics of the photon number and 
time distributions follows, essentially, a Poissonian law that severely restricts the range 
of practical applications of entangled photon sources based on PDC, e.g., for some 
quantum cryptography protocols [5]. 

In this paper, we propose a cavity Quantum Electrodynamics [HI [71 El [9l [101 dH 
[T2] (cavity-QED) implementation to deterministically generate polarization-entangled 
photons pairs that presents three important features from a practical point of view: (i) 
it operates with high fidelity even in the bad cavity limit; (ii) it is very robust under 
fluctuations of the system parameters since it is based on adiabatically following an 
energy eigenstate; and (iii) the initial field state of the system is the simplest in optical 
cavity-QED, namely the vacuum state for all cavity modes. 

Recently there have been several cavity-QED proposals for generating entangled 
photon pairs coupling a single atom to two e.m. modes of a single optical cavity or 
even by means of two optical cavities jlSl jlH [151 [IB]. Common to all these cavity- 
QED implementations is the requirement to operate in the strong coupling regime while 
the work that we present here is, to our knowledge, the first cavity-QED proposal for 
generating entangled photon pairs that operates even in the bad cavity limit, i.e., when 
the lifetime of the photon in the cavity, given by the transmission rate of the cavity 
mirrors, is much smaller that the typical time that the vaccum modes of the cavity need 
to produce a single quantum Rabi oscillation. Note also that a single-atom source in 
a cavity-QED setup has been considered for engineering entangled many-photon pulses 
consisting of a sequence of non-overlapping one-photon wave packets [IT] . Notably, this 
sequential generation has been completely characterized showing that the attainable 
states correspond to the hierarchy of the so-called matrix-product states [18] . 

The cavity-QED proposal that we discuss here to deterministically generated 
polarization entangled photon pairs is based on adiabatically following an energy 
eigenstate of the complete system atom-cavity modes that initially corresponds to an 
excited atom plus to empty cavity modes and, eventually, to the atom in its internal 
ground state and the two-photons entangled in their polarization degree of freedom. 
In Section 2, we introduce the physical scheme under investigation and derive the 
corresponding Hamiltonian. The coherent dynamics yielding the entangling protocol 
is discussed in Section 3. In section 4, we investigate the role of decoherence, i.e., 
spontaneous atomic decay and photon detection, by means of the Monte Carlo Wave 
Function approach. Section 5 is devoted to the characterization of the cavity-QED 
source entanglement capability while some practical considerations are briefly discussed 
in Section 6. Finally, in Section 7 we summarize the main results of the paper. 
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Figure 1. (a) Four-level atomic system: \Fmp) denotes the atomic levels, gia (i = 1, 2) 
(a = ±) account for the vacuum Rabi frequencies of each polarization mode, and Ajq, 
(a = ±) are the cavity detunings from the corresponding atomic transitions, (b) 
Manifold of atom- field states involving the initial state |/) = |0o) (8 |ri). (c) The same 
manifold as in (b) in the interaction picture and under the Raman resonance condition 
Ai+ = Ai_ and A2+ = A2-, and gi+ — gi-{— gi). For the definition of \B), \D), and 
\E^) see Eqs. (gl)-®- 



2. Physical framework 

The system under investigation, sketched in Fig. 1(a), is composed of two longitudinal 
cavity modes uoi and 002-, presenting polarization degeneracy, and two atomic transitions 
F = Q ^ F' = 1 and F' = 1' F" = 0" in a ladder configuration. As usual, F 
represents total angular momentum and mp its projection along the quantization axis. 
In what follows we will use the notation \Fmp) for the atomic state. The Hamiltonian 
of the system is composed of the free Hamiltonian of the atom and the e.m. cavity 
modes -ffatom and Hcav, respectively, and the interaction Hamiltonian in the rotating 
wave approximation Hi {h = 1), i.e.,: 

Ht = -f^atom + -f^cav + Hi (l) 
H^tora = ^2- |1-1> (1-1 1 + ^2+ |l'l) (I'll + (^1+ + UJ^-) \Oo) (OqI (2) 

-f^cav = ^ ^ i^ia (al^aia^ (3) 

i=l,2 0=+,- 

= XI XI 3ia{t) (^al^Sia + a^aSl^^ , (4) 
i=l,2 a=+,- 

where i = 1,2 denotes the two longitudinal modes and a = ± referes to the two circular 
orthogonal polarizations. The energies of the atomic states are given as a function of 
the electric dipole transition frequencies u!i± {i = 1, 2), with state |0q) being the zero of 
energies. gia{t) is the corresponding time-dependent vacuum Rabi frequency for each 
polarization mode. (ai±) is the photon creation (annihilation) operator for each 
mode, and 5i+ = |l'_i)(Oo|, 5i_ = |l'i)(Oo|, ^2+ = |0'o')(i;|, S2- = |0'o')(l'_il are atomic 
lowering operators. Detunings are defined as Ai± = ui — uji± and A2± = UJ2 — uj2±- 
For the sake of simplicity, we will assume equal vacuum Rabi frequencies for the two 
polarization states of each longitudinal cavity mode , i.e., gij^ = gi- = gi. 
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3. Entangling mechanism 

We take the initial state of the system to be = 0)) = |Oo)®|f2) with = |r2i)(8)|r22), 
being \ the vacuum state of mode i = 1,2. The coherent evolution of the system will 
remain in the space spanned by the five states of the manifold shown in Fig. 1(b). 
This evolution is, in general, very involved except for the Raman resonant case, 
Ai+ = Ai_(= Ai), and A2+ = A2-(= A2) where the interactions can be reduced 
to that of a three-level system. In order to show this, we will consider the following 
basis: 

|/) = |Oo) ® (5) 

V2\B) = (^S^+a\+ + Sl-a\-)\I), (6) 

V2\D) ^ (5i+ai+ - 5i_al_) |/), (7) 

V2\E^) = (^S2-al_Si+a\^ ± 52+4+5i_al_) |/). (8) 

\B) and \ D) are the so-called bright and dark state combinations of the atomic states 
|li) and |l-i) and the two circularly polarized states of mode Ui. \E^) correspond to the 
Bell states of two orthogonally polarized photons, one in each longitudinal mode, with 
the atomic state |0q) factorising. Under the Raman resonance condition mentioned 
above, it is straightforward to check that the interaction picture Hamiltonian of the 
system, H^, has the following off-diagonal matrix elements: 

{D\H{\I) = {DlHilE-') = {B\Hi\E-) = 0, (9) 
{B\Hi\I) =V2g,e~^^^\ (10) 
{BmE+) = {D\H{\E-)=g,e^^-'. (11) 

The coupling chain thus reduces to that of the three- level system shown in Fig. 1(c): 

|/) ^ \B) ^ \E+). 

Under the two-photon resonance condition Ai = — A2, one of the energy eigenstates 
of this three-level system is 

|A(^)) =cos^|/) -sin^|E+), (12) 

with tan^^(t) = \/2gi{t) / g2{t). For 6* = 0, Eq. (IT^ corresponds to the initial state of 
the system: |A(0)) = |/). Following the STIRAP technique [201 |2ll |22] , the adiabatic 
change of 6 from {gi = 0, (72 7^ 0) to 7r/2 {gi 7^ 0, (72 = 0) transfers the population from 
|/) to the entangled state \E^) without populating the intermediate state \B). Since 
|/) contains no cavity photons, the leakage of two photons through the cavity mirrors 
indicates the success in generating state \E^) which, in turns, shows the possibility of 
operating with high fidelity even in the bad cavity limit, i.e., for n > g. The STIRAP 
technique requires the counterintuitive interaction of the four level atom with the cavity 
mode UJ2 and later on, and with an appropriate temporal overlapping, with mode ui. 
Notably, this scheme is very robust under fluctuations of the experimental parameters 
provided the adiabaticity condition is satisfied during all the interaction process. 
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4. Monte Carlo Wave Function simulations 

To numerically investigate the previous proposal in the presence of dissipation, we 
will use next the Monte Carlo Wave Function (MCWF) formalism [23j. In this 
approach, the time evolution of the wave fucntion of a single quantum system, a so- 
called quantum trajectory, is calculated by integrated the time-dependent Schr'odinger 
equation using an effective non-Hermitian Hamiltonian. Incoherent process, such as 
spontaneous emission or photon detection, are incorporated as quantum-jumps ocurring 
at random times. Thus, a quantum trajectory consists of a series of coherent evolution 
periods separated by quantum-jumps. Notably, the MCWF formalism is equivalent to 
the density-matrix formalism but provides better insights into the underlying physical 
mechanisms [Ml [25l [26] and accordingly we will use it here. For the system under 
investigation, the non-hermitian Hamiltonian is 

He{i = H[-Y^ (i^Sj^Sia + i'^al^aiA, (13) 

where k accounts for the cavity-photon emission rate and the eventual photon-detection, 
and r corresponds to the spontaneous atomic decay rate. For simplicity, we have taken 
equal cavity and atomic decay rates for all cavity modes and for the four atomic 
transitions, respectively. As usual in the MCWF formalism, to calculate a single 
quantum trajectory a new pseudorandom number e is used at each interval time dt 
to decide whether a quantum jump occurs. In case it doesn't occur, the wavefunction 
has to be renormalized to ensure unitary evolution. 

After averaging over many quantum trajectories, we have obtained the probabilities 
associated to all possible photon-emission events shown in Figs. [2|^a) and [2|^c) for 
K = O.lg and F = 0.01^, and Figs. Mh) and EJ^d) for k = g and F = 0.01^. In 
the four cases, we have considered properly overlapped temporal Gaussian profiles for 
the interaction strengths. 

As it is inferred from Fig. EJ^b) and[2t^d), even in the bad cavity limit, the probability 
of generating an entangled pair of photons (corresponding to column (i)), is larger than 
the probability of creating a separable state of two photons emitted from the cavity 
(column (ii)). The later process is due to the photon-cavity decay from state \B), 
i.e., it originates from diabatic processes. By speeding up the sequence. Figs. |2](c- 
d), the population of \B) and the probability of generating two non-entangled photons 
will increase, while the probability of spontaneous emission processes, corresponding to 
columns (iii) and (iv), will reduce. In general, the time duration of the sequence will be 
optimized to maximize the probability given in colum (i) but also the ratio between the 
probabilities (i) and (ii). 

The previous results can be largely improved by operating out from the single- 
photon resonance condition while maintining the two-photon resonance one, i.e., Ai = 
— A2 7^ 0. Thus, Fig. 3 shows the results averaging over many MCWF simulations 
for the following set of parameters k = 2g and F = O.Olg, yielding (a) F = 0.79 



Generation of entangled photon pairs in the bad cavity limit 6 

1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

0.0 
1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

Figure 2. Probabilities for (i) generating a maximally entangled photon pair; (ii) 
generating two non-entangled cavity photons; (iii) emitting one cavity photon and one 
spontaneously emitted photon; (iv) emitting two photons by means of spontaneous 
emission. Parameters are k = O.lg and F = O.Olgf for (a) and (c), and k = g 
and F = O.Ol.g for (b) and (d). In all cases Ai = A2 = 0. Results arc obtained 
averaging over many Monte Carlo wave function simulations. The time duration of 
the entanglement procedure is gt= 100 for (a) and (b) and gt = 30 for (c) and (d). 
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Figure 3. Probabilities for events (i)-(iv). Parameters are Ai = — A2 = 5g for (a), 
Ai = — A2 = IQg for (b) and Ai = — A2 = 15g for (c). In all cases k = 2g and 
F = O.Olg and the time duration of the entanglement procedure is gt = 100. Results 
are obtained averaging over many Monte Carlo wave function simulations. 

for Ai = -A2 = 5g, (b) F = 0.96) for Ai = -A2 = lOg and (c) F = 0.99 for 
Ai = -A2 = I5g. 

5. Chciracterization of the source 

By means of coincidence photodetection, one can keep only those events for which one 
photon is detected from each cavity mode and evaluate the fidelity F — {E~^\p\E~^) where 
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Figure 4. The CHSH parameter S and the FideUty as a function of the cavity 
decay rate. The parameters are F = O.Olg, gt = 100 and Ai — — A2 = 7g for the 
pointed Une and Ai = — A2 = Ibg for the sohd line. 
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Figure 5. (a) Entanglement capability 5 as a function of the deviation from the two- 
photon resonace condition and the delay time gtd between the time-dependent coupling 
strengths of modes 1 and 2. Parameters are n = g^V — O.Olg, Ai+ = Ai_ — Ai, 
A2+ = A2- = A2 and the width (FWHM) of the coupling pulses is gt^ = 27. The 
time duration of the entanglement procedure is, gt = 100. 

p is the corresponding reduced density matrix. In addition, to quantify the entanglement 
capabihty of the source we have used the S parameter of the CHSH inequahty [27] 
{S = \/2 for any separable state, S = 2\/2 for maximally entangled states, and S < 2 
for local hidden variable theories). Thus, for the simulations corresponding to Fig. 3 
one obtains (a) (F, S) = (0.79, 2.25) for Ai = -A2 = 5g, (b) (F, S) = (0.96, 2.73) for 
Ai = -A2 = 10^ and (c) {F, S) = (0.99,2.81) for Ai = -A2 = 15g. For the set of 
parameters F = O.Olg and gt = 100, Fig. 4 shows the fidelity F and the entanglement 
capibility 5* for Ai = — A2 = 7g (dotted line) and Ai = — A2 = 15g (sohd line), as 
a funcion of the cavity losses k. It is clearly seen that the cavity-QED proposal here 
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discussed should operate with high fidehty well inside the bad cavity regime, e.g., for 
K ~ 3g, provided the the two cavity modes are far from single photon resonance but 
fulfill the two-photon resonance condition. 

In Fig. ([5]) the entanglement capability S is shown as a function of the deviation 
from the two photon resonance condition and the delay time between the coupling 
strengths gi and g2- Note that as soon as the two photon resonance condition is broken 
the STIRAP procedure fails and the S parameter decreases. On the other hand, the S 
parameter exhibits a robust behavior against the variation of the delay time between 
the two interactions as it is seen in the horizontal axis of Fig. ([5]). 

6. Some practical considerations 

A suitable atomic element with the configuration F = ^ F' = 1 and F' = 1' ^ 
F" = 0" needed in our proposal, is calcium with its cascade J = 0-^J=1^J = 
transition V ^'^o ^ 45^^^^! ^ ^Sq at A = 551.3 nm and A = 422.7 nm, respectively 
[281 [29] . Alternative one could also use aJ = l^J=l^J = configuration 
like the 7^5*1 —>■ 6^ Pi Q^Sq of the ^'^''Hg [30]. In both cases, one needs to perform 
first a double excitation process such that the atom enters the cavity-QED setup in 
its excited state. Note that as the initial state can decay by spontaneous emission, if 
one pretends to obtain a close to 100% success probability in generating the entangled 
photon pair, then the spontaneous population decay rate should be much smaller than 
the inverse of the total time needed for the STIRAP process. In fact, our simulations 
(see Figs. 2 to 5) include not only cavity decay of photons through the mirrors but 
also the spontaneous atomic decay from the optical transitions. For the simulations of 
these figures, the spontaneous atomic decay rate has been taken smaller but close to the 
inverse of the total time needed for the STIRAP process. However, if the spontaneous 
emission rate becomes on the order or much larger than the total time for the STIRAP 
process, then the success probability will obviously decrease. If so, one could perform 
postselection of events, i.e, two photon coincidence measurements, and still achieve a 
fidelity close to 100% in generating entangled photon pairs. 

Finally, note that the temporal control of the interactions, i.e., the switching 
on/off of the interaction of the atomic transition with the corresponding mode, can 
be implemented by detuning the atomic transitions from he corresponding modes via a 
uniform electric (or magnetic) field yielding the corresponding Stark (or Zeeman) shift, 
or even by the atomic light-shift induced by an external laser. 

7. Conclusion 

To sum up, we have discussed a robust and efficient cavity-QED proposal for 
the deterministic generation of polarization-entangled photon pairs. The complete 
entanglement procedure is based on a two mode STIRAP process which allows the 
cavity-QED source to operate with high fidelity even in the bad cavity limit. Even in 
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this regime, the obtained fidelities yield quantum correlations between photon pairs well 
above from those that can be obtained from local hidden variable theories. 
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